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1 Lyapunov coefficients 



The beginning of this section is a review of the method found in [3] and in [I] 
for the calculation of the first and second Lyapunov coefficients. The calcu- 
lation of the third Lyapunov coefficient can be found in [6]. The calculation 
of the fourth Lyapunov coefficient has not been found by the authors in the 
current literature. The extensive calculations and the long expressions for 
these coefficients have been obtained with the software MATHEMATICA 5 
[TT] and the main computational steps have been posted in the site |10j . 
Consider the differential equations 

x' = /(x,/i), (1) 

where x 6 W a and fj, e K m are respectively vectors representing phase vari- 
ables and control parameters. Assume that / is of class C°° in ¥L n x W 71 . 
Suppose (pQ) has an equilibrium point x = x at [i = /i and, denoting the 
variable x — x also by x, write 

F(x) = /(x,//o) (2) 

as 

F(x) = Ax + - B(x, x) + - C(x, x, x) + — D(x, x, x, x) + E(x, x, x, x, x) + 

^ K(x, x, x, x, x, x) + L(x, x, x, x, x, x, x) + M(x, x, x, x, x, x, x, x) (3) 

H N(x, x, x, x, x, x, x, x, x) + 0(||x| 1 10 ), 
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where A = / x (0, /x ) an d 



Bi(x,y) = £ 



£j 2/fc, (4) 

5=0 



Oi(x,y,z) = 2^ 



Xj Vk Zh (5) 
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(x,y,z ' u) ^£w^ 



Xj Vk Zl U r , 



f=0 



#i(x,y,z,u,v) = ^ 



■ t ; 1 dfj d£fc <9£ r di p 

j,k,l,r,p=l J r 



Xj Uk Zl U r V p , 



(6) 



(7) 



e=o 



Ki(x,y,z,u, v, w) = ^2 



n 



Li(x,y,z,u, v,w,t) = ^ 

j /' I 



-J =l d£j d£ k d& d£ r d£ p d£ q 

d 7 F(o 



Xj y k Zl U r Vp W q , (8) 



d£j d£k d£i d£ r d£ p d£ q d£ h 



Xj y k ZlU r VpW q t h , (9) 



Mj(x,y,z,u, v,w,t,r) = ^ 



iVi(x,y,z,u,v,w,t,r,s) = > — — 

,3T : ^ ■ ■ ■ % 



y fc ^ U r Vp w q t h r a , (10) 



Xj Vk Zl U r Vp Wq t h T a S b , (11) 



e=o 



for % — 1, . . . , n. 



Suppose (xo, Ho) is an equilibrium point of ([I]) where the Jacobian matrix 
A has a pair of purely imaginary eigenvalues A2,3 = ±io>o, u > 0, and admits 
no other eigenvalue with zero real part. Let T c be the generalized eigenspace 
of A corresponding to A2,3- By this is meant that it is the largest subspace 
invariant by A on which the eigenvalues are A2,3- 

Let p, q G C n be vectors such that 



Aq = iu q, A T p = -iu p, (p, q) = ^Pi ft = 1, 



(12) 



i=\ 



where A T is the transposed of the matrix A. Any vector y G T c can be 
represented as y = wq + wq, where w = (p, y) G C. The two dimensional 
center manifold can be parameterized by w, w, by means of an immersion of 
the form x = H(w,w), where H : C 2 — > 1R™ has a Taylor expansion of the 
form 



H(w,w) = wq + wq + — ^— h jk w j w h + 0(1 

v.k\ 

2<J+fc<9 J 



w 



10 \ 



(13) 
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with hj k G C n and hj k = h k j. Substituting this expression into (JTJ) we obtain 
the following differential equation 

H w w' + = F(H(w,w)), (14) 

where F is given by (j2J). 

The complex vectors h%j are obtained solving the system of linear equa- 
tions defined by the coefficients of (fT4"j) . taking into account the coefficients 
of F, so that system (fT4"j) . on the chart w for a central manifold, writes as 
follows 

1 1 11 

u;' = iuj w+- G 2 iw\w\ +— G 32 w\w\ 4 +— G 43 w\w\ 6 +— — - G 54 w\w\ 8 +0(\w\ 10 ), 

Z LZ 144 ZooU 

with Gjk € C. 

Solving for the vectors /iy the system of linear equations defined by the 
coefficients of the quadratic terms of (|14p . taking into account the coefficients 
of F in the expressions (jSJ) and fll]), one has 

ftn^-A-^^g), (15) 

/ i20 = (2^ / n -A)- 1 J B( ? ,g), (16) 

where J n is the unit n x n matrix. Pursuing the calculation to cubic terms, 
from the coefficients of the terms w 3 in (fl~4"l) follows that 

h 30 = (3iu Q I n - A)' 1 [3B(q, h 20 ) + C(q, q, q)\ . (17) 

From the coefficients of the terms w 2 w in ([14"]) one obtains a singular 
system for h 2 \ 

(iu I n - A)h 21 = C(q, q, q) + B(q, h 20 ) + 2B(q, h n ) - G 21 q, (18) 
which has a solution if and only if 

(p, C(q, q, q) + B(q, h 20 ) + 2B(q, h u ) - G 2l q) = 0. 
The first Lyapunov coefficient l\ is defined by 

h = ~ReG 21 , (19) 
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where 



G 2 i = (jp, H21) , and H 2 i = C (q, q, q) + B(q, h 20 ) + 2B(q, h 



11; 





The complex vector h 2 \ can be found by solving the nonsingular (n + 1)- 
dimensional system 



(20) 



with the condition (p, h 2 \) = 0. The procedure above can be adapted in 
connection with the determination of h 32 and 

For the sake of completeness, in Remark 11.11 we prove that the system 
(!20|) is nonsingular and that if (v, s) is a solution of (1201 with the condition 
(p, v) = then v is a solution of ( TT8l) . See Remark 3.1 of [6]. 

Remark 1.1 Write M. n = T c © T su , where T c and T su are invariant by A. 
It can be proved that y G T su if and only if {p, y) = 0. Define 

a = C(q, q, q) + B(q, h 20 ) + 2B(q, hu) - G 21 q. 

Let (v , s) be a solution of the homogeneous equation obtained from [2^1 . 
Equivalently 

(iu I n - A)v + sq = 0, (p, v) = 0. (21) 

From the second equation of l[21\) . it follows that v G T su , and thus (iuJo^n ~ 
A)v G T su . Therefore (p, {iuioI n — A)v) = 0. Taking the inner product of 
p with the first equation of (21\) one has (p, (iuJoI n — A)v + sq) = 0, which 
can be written as (p,(iuoI n — A)v) + s(p,q) = 0. Since (p,q) = 1 and 
(p, (iu I n — A)v) =0 it follows that s = 0. Substituting s = into the first 
equation of l[21\) one has (iuJ I n — A)v = 0. This implies that 

v = aq, aE C. (22) 
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But = (p,v) = (p,aq) = a(p,q) = a. Substituting a = into |HP it 
follows that v = 0. Therefore (v, s) = (0, 0). 
Let (v, s) be a solution of Equivalently 

(iu; I n — A)v + sq = a, (p, v) = 0. (23) 

From the second equation of (d3j) ; it follows that v G T su and thus (iu)oI n — 
A)v G T su . Therefore (p, (iu>oL n — A)v) = 0. Taking the inner product of p 
with the first equation of ( E^ ) one has (p, (iu> I n — A)v + sq) = (p, a) , which 
can be written as 

(p, (iu I n ~ A)v) + s(p, q) = (p, a). 

As (p, a) = 0, (p,q) = 1 and (p, {iuj^ln — A)v) = it follows that s = 0. 
Substituting s = into the first equation of [23^1 results (iu)oI n — A)v = a. 
Therefore v is a solution of l[T8\) . 

From the coefficients of the terms w 4 , w 3 w and w 2 w 2 in (fT4l) . one has 
respectively 

h m = {Aiu^In - A)~ l [35(^20, h 2 Q ) + 4B(g, h 30 ) + 6C(g, q, h 20 ) + -D(g, q, q, q)}, 

(24) 

/i 3 i = (2iu I n - A)- 1 [35(g, /i 21 ) + B(q, h 30 ) + 3B(h 20 , h u ) 

+3C(g, q, h u ) + 3C(q, q, h 20 ) + D(q, q, q, q) - 3G 21 h 20 ], (25) 

h 22 = -A~ l [D(q, q, q, q) + 4C(g, q, h n ) + C(q, q, h 20 ) + C(q, q, h 20 ) 

+2B{h u , /in) + 2B(q, h 21 ) + 2B(q, h 2l ) + B{h 2 o, h 20 )}, (26) 

where the term — 2hn{G 2 i + G 2 i) has been omitted in the last equation, since 
G 21 + G 21 = as h = 0. 
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Defining H32 as 

H 32 = 6B(h u ,h 21 ) + B(h 20 , h 30 ) + 3B(h 21 ,h 20 ) + 3B(q, h 22 ) 
+2B(q, h 31 ) + 6C(q, h n , h n ) + 3C{q, h 20 , h 20 ) + 3C{q, q, h 21 ) 
+6C(q, q, h 21 ) + 6C(q, h 20 , hn) + C(q, q, h 30 ) + D(q, q, q, h 20 ) 
+QD(q, q, q, h n ) + 3D(q, q, q, h 20 ) + E(q, q, q, q, q) 

—QG 2 ih 2 i — 3G21/121, 



and from the coefficients of the terms w 3 w 2 in (TH1) . one has a singular system 
for h 32 

(iu I n - A)h 32 = H 32 - G 32 q, (27) 
which has solution if and only if 

(p,H 32 -G 32 q) = 0. (28) 



where the terms —QG 2 \h 2 i — 3G 2 \h 2 i in the last line of (1271) actually does not 
enter in last equation, since (p, h 2 i) = 0. 

The second Lyapunov coefficient is defined by 

Z 2 = ^ReG 32 , (29) 

where, from ®, G 32 = {p,H 32 ). 

The complex vector h 32 can be found solving the nonsingular (n + 1)- 
dimensional system 



iu I n - A q\ / h 32 \ I H 32 - G 32 q 
P J 



(30) 



with the condition (p, h 32 ) = 0. 

From the coefficients of the terms w 4 w, w A w 2 and w 3 w 3 in f|T4l) . one has 
respectively 

h A1 = {3iu I n - A)- x [±B(hxi, h 30 ) + 6B(h 20 , h 2l ) + 4B(q, h 31 ) 



+B(q, h m ) + 12C(g, /in, h 20 ) + 6C(g, g, /i 21 ) + 4C(g, g, /i 30 ) (31) 
+3C(g, /i 20 , ^20) + 4D(g, q, q, h u ) + 6L>(g, q, q, h 20 ) 

+E(q,q,q,q,q) - 6G 21 h 30 ], 

h 42 = (2iu I n - A)- 1 [8 J B(/i 11 , h 31 ) + 6B(h 20 , h 22 ) + B(h 20 , h 40 ) 
+QB(h 21 ,h 21 ) + AB(h 2U h 30 ) + 45(g, /i 32 ) + 2B(q, h 41 ) 
+12C(h u , h u ,h 20 ) + 3C(h 20 , h 20 , h 2 o) + 24C(g, h u , h 21 ) 
+12C(g, /i 20 , M + 4C(g, /i 20 , M + 6C(g, g, /i 22 ) + 8C(g, g, /i 3 i) 

+8C(g, /i n , hw) + 12C(g, /i 20 , /i 2 i) + C(q, g, /i 40 ) (32) 
+ 12D(g, g, /in, /in) + 6D(q, g, /i 20 , /i 20 ) + 4D(g, g, g, /i 2i ) 
+12D(g, g, g, /i 21 ) + 24D(g, g, /in, /i 20 ) + 4D(g, g, g, /i 30 ) 
+3L>(g, g, h 20 , h 20 ) + E(q, g, g, g, /i 20 ) + 8£"(g, g, g, g, /in) 
+6£(g, g, g, g, /i 20 ) + K(q, q, q, q, g, g) 
-4(G 32 /i 20 + 3G 2 i/i 31 + G 21 h 31 )\, 

h 33 = -A-^Bihu, h 22 ) + 3B(h 20 , h 31 ) + 3S(/i 2 o, h 31 ) + 9B(h 21 ,h 21 ) 
+B(h 30 , h 30 ) + 3B(q, h 32 ) + 35(g, h 32 ) + 6C(/i n , /in, /in) 
+9C(/in, h 20 , h 20 ) + 18C(g, /in, h 21 ) + 3C(g, /i 20 , /i 30 ) 
+9C(q, h 20 , h 21 ) + 3C(g, g, /i 31 ) + 9C(g, g, h 22 ) + 18C(g, /i n , /i 2i ) 
+9C(g, /i 20 , /i 2i ) + 3C(g, h 20 , h 30 ) + 3C(g, g, /i 3i ) + 9D(g, g, /i 20 , hu) (33) 
+D(g, g, g, h 30 ) + 9D(g, g, g, h 21 ) + 18L>(g, g, /i n , h u ) 
+9D(q, g, /i 20 , /i 20 ) + 9D(g, g, g, /i 21 ) + 9£>(g, g, h u ,h 20 ) 
+SE(q, g, g, g, /i 20 ) + 9E(q, g, g, g, /in) + 3E(q, g, g, g, h 20 ) 
+K(q, q, g, g, g, g) - 3(G 32 + G 32 )h n - 9(G 21 + G 21 )/i 22 ]. 

Defining H 43 as 

W43 = l2B(h u ,h 32 ) + 6£(/i 20 , M + 3B(h 20 , h 41 ) 
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+ 18B(h 21 ,h 22 ) + l2B{h 21 ,h 31 ) + 4B{h 30 , hi) + 5(^30, h4o) 
+4B{q, h 33 ) + W(q, h 42 ) + 36C(h u , h n , h 21 ) + 36C{h n ,h 20 , h 21 ) 
+12C(h u ,h 20 , h 30 ) + 3C(/i 20 , h 20 , h 30 ) + l8C(h 20 , h 20 , h 21 ) 
+mC(q, h n ,h 22 ) + 12C(g, h 20 , h 31 ) + 12C(q, h 20 , h 31 ) 
+36C(g, h 21 , h 21 ) + 4C(g, h 30 , h 30 ) + 6C(g, g, h 32 ) 
+12C(g, g, /i 32 ) + 24C(g, h n , h 3l ) + 18C(g, h 20 , h 22 ) 
+3C(g, ^. 2 o, /i 4 o) + 18C(g, /i 2 i, ^21) + 12C(g, h 2X , h 30 ) 
+3C(g, g, /i 4 i) + 24D(g, /in, /in, h u ) + 36D(g, /in, /i 20 , ^20) 
+36L>(g, g, hn,h 2 i) + 6-D(g, g, /i 2 o, ^30) + 18£>(g, g, /i 20 , ^21) 
+4D(g, g, g, h 3 i) + l8D(q, q, g, h 22 ) + 72D(q, g, hu,h 2 i) 
+36L>(g, g, /i 20 , h 2 l) + 12£>(g, g, ^20, h 30 ) + 12£>(g, g, g, h 31 ) 
+36D(g, /in, ^11, ^20) + 9£>(g, /i 20 , ^20, ^20) + 12£>(g, g, /in, /130) 
+ 18Z)(g, g, h 20 , h 21 ) + £>(g, g, g, /i 40 ) + 12E'(g, g, g, h u ,h 20 ) 
+E{q, g, g, g, /130) + 12£(g, g, g, g, ^21) + 36£(g, g, g, /in, /in) 
+18£'(g, g, g, h 20 , h 20 ) + l8E(q, g, g, g, h 21 ) + 36E(q, g, g, h u ,h 20 ) 
+4£(g, g, g, g, h 30 ) + 3E(q, g, g, h 20 , h 20 ) + 3£T(g, g, g, g, g, /i 20 ) 
+ 12K(g, g, g, g, g, /in) + 6K(q, g, g, g, g, h 20 ) + L(g, g, g, g, g, g, g) 
-6{2G 32 h 21 + G 32 /i 2 i + 3G 2 i/i 32 + 2G 21 h 32 ), 

and from the coefficients of the terms w 4 w 3 in ffT4l) . one has a singular system 
for /i 43 

(iuj In- A)h 43 =H A3 - G 43 q (34) 

which has solution if and only if 

(p,H i3 -G 43 q) = 0. (35) 

where the terms — 6(2G 32 /i 2 i + G 32 h 2 i + 3G 2 i^3 2 + 2G 2 i/i 32 ) appearing in the 
last line of equation ( 1341) actually do not enter in the last equation, since 
(p, h 21 ) = and (p, h 32 ) = 0. 

The third Lyapunov coefficient is defined by 

/ 3 = -I-ReG 43 , (36) 
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where, from (13"5|) . G43 = (p, W43). 

The complex vector ^,43 can be found solving the nonsingular (n + 1)- 
dimensional system 



(37) 

p 

with the condition (p, h 43 ) = 0. 

Defining 7-^54 by the expression below 

2ftB(ftu, /143) + lOB(hm, hs) + QB(ho, h 52 ) + 40B{h 21 ,h 33 ) + 30B(hi 
60B(h 22 , h 32 ) + 10B{h 30 , hi) + 4B(h 30 , hsi) + 4QB(h 31 ,hi) + 20B(h 3 i 
5B(h 40 , hi) + B(h 40 ,h 50 ) + 5B(q, liu) + 4B(q, h 53 ) + 120C(h n ,h n 



A q\ 


II 


( h 43 \ 




0) 




K s J 


H 




60C(fc u , ho, M + 360C(/i u , h 21 ,h 22 ) + 240C(hn,hi,h 31 ) + 8QC(hu 
20C(h u , ho, M + l20C{h 20 , h xx ,hi) + l5C{h 20 , h 2Q , hi) + 60C(h 20 
!20C(h 20 , h 21 ,h 3 i) + lSOC(h 20 , hi,h 22 ) + WC(h 20 , h 3Q , ho) + 40C(h 20 
3C(ho, ho, h 50 ) + 120C(ho, h 21 ,h 31 ) + 30C(ho, hi, h 40 ) + 60C(ho 
180C(h 21 , h 2 i, h 2 i) + 60C(/i2i, h 2 i, Ii 3 q) + 40C(h 3 o, h 2 \, h 3 o ) + 8QC{q 
30C(q, h 2Q , hi) + 30C7(q, ho, h 42 ) + 120C(q, h 21 , h 3 l) + l2QC(q 
90C(q, hn, h 22 ) + 20C(q, h 30 , hi) + 20C(q, h 3 0, M + &0C(q, h 31 ,h 3 l) + 5C(q 
10C(q, q, hz) + 20C(q, q, h 43 ) + 60C(q, h 11 ,h 42 ) + 40C(g, h 20 , h 33 ) + \2C{q 
120C(g, h 21 ,h 32 ) + 60C(g, hi,h 41 ) + 40C(q, h 30 , hi) + 4C(q 
120C(g, h 3 i, h 22 ) + 20C(q, h 4Q , hi) + 6C($, q, h 52 ) + 240D(h u , h n 
120D(h n , h n ,ho, h 30 ) + 360D(/i 20 , hn, hn, hi) + 360D(h 20 , hn 
60D(h 20 , h 2 o, hn,ho) + 90D(h 20 , h 20 , ho, hi) + 30D(h 20 , ho 
36(XD(g, h n ,hn,h 22 ) + 240D(q, hn, ho, h 3 i) + 720D(q, hn 
80D(q, hn,h 30 , ho) + 24QD(q, h 2Q , hn,hi) + VaD{q, h 20 
180D(g, hio, ho,h 22 ) + 120D(q, h 20 , h 21 ,ho) + 180D(g, h 20 
15D(q, ho, ho,h 4Q ) + 180D(q, ho, ho, hi) + 120D(g, ho 
l20D(q, q, hiuhl) + 30D(q, q, h 20 , hi) + 60D(q, q, ho, h 32 ) + 1201% q 
1802% q, hi,h 22 ) + 102%, q, h 30 , ho) + 402%, q, ho, h 31 ) + 10D(q 
402%, q, q, h 33 ) + 240D(q, q, h n , h 32 ) + 120D(q, q, h 20 , hi) + 60£>(g, q 
360D(q, q, h 21 ,h 22 ) + 240D(q, q, hi, h 31 ) + 80D(q, q, h 30 , hi) + 20D{q, q 
'S0D(q, q, q, h 42 ) + 240D(q, hn,h n , h 31 ) + 60D(q, hn, ho, h ) + 360D(q, h u 
240D(q, h n ,h 30 , hi) + 360D(q, h 20 , h n ,h 22 ) + 60D(q, h 20 
120D(q, h 20 , h 20 , h 31 ) + 36QD(q, h 20 , h 21 ,h 21 ) + 40D(q, h 20 
120D(q, ho, h 30 , h 21 ) + 60D(q, q, hn,h 41 ) + 60D(g, q, h 2Q ,h 32 ) + 6D(q, q 



h 3 o 
ho 
ho 
h 3 o 
hn 
hi 
h 4 o 
ho 
ho 
hn 
ho 
ho 
h 2 i 
h 20 
hi 
h 3 o 
h 2 i 

q,q 

ho 
ho 
hi 
h 20 
ho 
h 20 
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120£>(<7, q, h 21 ,h 31 ) + 30D(q, q, hi,h4o) + GOD(q, q, h 30 , h 22 ) + 4D(q, q, q, h 51 ) + 
120E(q, hn, hn.hu, h n ) + 360E(q, h 20 , hn,h n ,h 20 ) + 45£%, h 2 o, h 20 , h 20 , h 20 ) + 
360E(q,q, hn, hn, h 2 \) + 360E(q, q, hn, h 2 o, h 2 i) + 120E(q, q, h 2 o, hn, h 3 o) + 
180E(q, q, h 20 , h 20 , h 21 ) + 30E(q, q, h 20 , h 20 , h 3Q ) + 80E(q, q, q, hn, h 3 i) + 
10E(q, q, q, h 20 , h ia ) + 60E(q, q, q, h 20 , h 22 ) + 40£%, q, q, h 21 ,h 3a ) + 60E(q, q, q, h 21 , h 21 ) + 
5E(q, q, q, q, h 41 ) + AOE(q, q, q, q, h 32 ) + 360E(q, q, q, hn, h 22 ) + 120 E(q, q, q, h 20 , h 31 ) + 
120 E(q, q, q, h 20 , h 3 i) + 360 E(q, q, q, h 21 , h 21 ) + AOE{q, q, q, h 30 , h 30 ) + 60E(q, q, q, q, h 32 ) + 
720E(q, q, hn, hn,h 21 ) + 2A0E(q, q, hn, h 20 , h 30 ) + 720E(q, q, h 20 , hn, h®i) + 
60E(q, q, h 2Q , h 20 , h 30 ) + 360 E(q, q, h 20 ,h 2a , h 21 ) + 240E(q, q, q, hn, h 31 ) + 
180E(q, q, q, h 20 , h 22 ) + 30E{q, q, q, h 20 , h i0 ) + 180E(q, q, q, h 21 , h 21 ) + 120E(q, q, q, h 30 , h 21 ) + 
20E(q, q, q, q, h A1 ) + 240E(q, h 20 , hn, hn,h n ) + 180 E(q, h 20 , h 20 , hn, h 2 o) + 
120E(q, q, hn, hn, h 30 ) + 360E(q, q, h 20 , hn,h 2 %) + 90E(q, q, h 20 , h 20 , h 21 ) + 
60E(q, q, h 20 , h 20 , h 30 ) + 20E(q, q, q, hn, /i4o) + 40-E l (g, q, q, h 20 , h 31 ) + 40E(q, q, q, h 3Q , h 21 ) + 
E(q, q, q, q, h 50 ) + 120K(q,q,q,hn,hn,h 20 ) + 30K{q,q,q,h 2Q ,h 2a , h 20 ) + 
20K(q, q, q, q, hn, h 30 ) + 30K(q, q, q, q, h 20 , h 21 ) + K(q, q, q, q, q, h 40 ) + 20K(q, q, q, q, q, h 31 ) + 
2A0K(q,q,q,q,hn,h 21 ) + A0K(q,q,q,q, h 20 ,h 30 ) + 120K(q,q,q,q, h 20 ,h 21 ) + 
60K(q, q, q, q, q, h 22 ) + 24QK{q, q, q, hn,hn, h n ) + 360K (q, q, q, h 20 , hn, h 20 ) + 
360K(q, q, q, q, hn, h 21 ) + 180K(q, q, q, q, h 20 , h 21 ) + 60K(q, q, q, q, h 20 , h 30 ) + 
AOK{q, q, q, q, q, h 31 ) + 360K(q, q, q, h 2Q ,hn, hn) + 90K(q, q, q, h 20 , h 20 , h 20 ) + 
80K(q, q, q, q, hn,h 3Q ) + 120K(q, q, q, q, h 20 , h 21 ) + 5K(q, q, q, q, q, hio) + 
60K(q, q, q, h 20 , h 2Q , hn) + WK(q, q, q, q, h 2Q , h 30 ) + 3L(q, q, q, q, q, h 20 , h 20 ) + 
4£(?, q, q, q, q, q, ho) + Q0L(q, q, q, q, q, hn,h w ) + 30L(g, q, q, q, q, q, h 21 ) + 
120L(q, q, q, q, q, hn,h n ) + Q0L(q, q, q, q, q, h 20 , h 2Q ) + A0L(q, q, q, q, q, q, h 21 ) + 
120L(q, q, q, q, q, h 20 , h n ) + 10L(q, q, q, q, q, q, h 3a ) + 15L(q, q, q, q, q, h 20 , h 20 ) + 
6M(q, q, q, q, q, q, q, h 20 ) + 20M(q, q, q, q, q, q, q, h n ) + 10M(q, q, q, q, q, q, q, h 20 ) + 

N{q,q,q,q,q,q, q, q, q), 

and from the coefficients of the terms w b w^ in (|14|) . one has a singular system 
for h 5i 

(iuj I n - A)h 54 = H 54 - G 54 q (38) 
which has solution if and only if 

(p,H 54 -G 54 q) = 0. (39) 
The fourth Lyapunov coefficient is defined by 

'*=2S) ReG «> (40) 
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where, from (I3"§j) . G 54 = (p,H 54: ). 

Remark 1.2 Other equivalent definitions and algorithmic procedures to write 
the expressions for the Lyapunov coefficients lj,j = 1,2,3,4, for two dimen- 
sional systems can be found in Andronov et al. [1] and Gasull et al. /![/, 
among others. These procedures apply also to the n-dimensional systems of 
this work, if properly restricted to the center manifold. The authors found, 
however, that the method outlined above, due to Kuznetsov |5j \4}j, requiring 
no explicit formal evaluation of the center manifold, is better adapted to the 
needs of long calculations in J3], 0, \7jj, and for that in J2J/, where n = 3. 

A H op f point (xo,/^o) is an equilibrium point of ([1]) where the Jacobian 
matrix A = / x (x ,/io) has a pair of purely imaginary eigenvalues A2,3 = 
±iuoo, > 0, and admits no other critical eigenvalues — i.e. located on the 
imaginary axis. At a Hopf point a two dimensional center manifold is well- 
defined, it is invariant under the flow generated by ([1]) and can be continued 
with arbitrary high class of differentiability to nearby parameter values. In 
fact, what is well defined is the oo-jet — or infinite Taylor series — of the 
center manifold, as well as that of its continuation, any two of them having 
contact in the arbitrary high order of their differentiability class. 

A Hopf point is called transversal if the parameter dependent complex 
eigenvalues cross the imaginary axis with non-zero derivative. In a neigh- 
borhood of a transversal Hopf point — HI point, for concision — with li ^ 
the dynamic behavior of the system ([T]), reduced to the family of parameter- 
dependent continuations of the center manifold, is orbitally topologically 
equivalent to the following complex normal form 

w' = (rj + ioo)w + li%u\w\ 2 , 

w G C, rj, uj and l\ are real functions having derivatives of arbitrary high 
order, which are continuations of 0, ujq and the first Lyapunov coefficient at 
the HI point. See (3J. As li < {l\ > 0) one family of stable (unstable) 
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periodic orbits can be found on this family of manifolds, shrinking to an 
equilibrium point at the HI point. 

A Hopf point of codimension 2 is a Hopf point where l\ vanishes. It is 
called transversal if n = and li = have transversal intersections, where 
7] = r](/j,) is the real part of the critical eigenvalues. In a neighborhood 
of a transversal Hopf point of codimension 2 — H2 point, for concision — 
with I2 7^ the dynamic behavior of the system ([1]), reduced to the family 
of parameter-dependent continuations of the center manifold, is orbitally 
topologically equivalent to 

w' = (r) + iu )w + tw\w | 2 + l2w\w\ 4 , 

where 77 and r are unfolding parameters. See [3]. The bifurcation diagrams 
for I2 7^ can be found in [3], p. 313, and in [9]. 

A Hopf point of codimension 3 is a Hopf point of codimension 2 where 

12 vanishes. A Hopf point of codimension 3 point is called transversal if 
rj — 0, l\ = and I2 = have transversal intersections. In a neighborhood 
of a transversal Hopf point of codimension 3 — H3 point, for concision — 
with (3^0 the dynamic behavior of the system (JTJ), reduced to the family 
of parameter-dependent continuations of the center manifold, is orbitally 
topologically equivalent to 

w' = (r) + iujo)w + tw\w\ 2 + uw\w\ A + l3w\w\ 6 , 

where 77, r and v are unfolding parameters. The bifurcation diagram for 

1 3 7^ can be found in Takens [9] and in [7]. 

A Hopf point of codimension ^ is a Hopf point of codimension 3 where 
£3 vanishes. A Hopf point of codimension 4 is called transversal if 77 = 0, 
l\ — 0, I2 — and I3 = have transversal intersections. In a neighborhood 
of a transversal Hopf point of codimension 4 — H4 point, for concision — 
with I4 / the dynamic behavior of the system (JTJ), reduced to the family 
of parameter-dependent continuations of the center manifold, is orbitally 
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topologically equivalent to 

w' = [j] + iuj Q )w + tw \w\ 2 + z/w|ty| 4 + crw\w\ 6 + Uw\w\ 8 , 
where r/, r, v and cr are unfolding parameters. 
Theorem 1.3 Suppose that the system 

x = /(x, /x), x = (x, y, z), fj, = a, ac, e) 
/ias t/ie equilibrium x = for (jl — wif/i eigenvalues 

A 2 ,3(» = »7(a*) ± 

where uj(0) = ujq > 0. For /i = i/ie following conditions hold 

r/(0) = 0, /i(0) =0, Z 2 (0) = 0, Z 3 (0) = 0, 

where li(fi), an d hit 1 ) are the first, second and third Lyapunov co- 

efficients, respectively. Assume that the following genericity conditions are 
satisfied 

1- ^(0) 7^ 0, where Z 4 (0) is the fourth Lyapunov coefficient; 
2. the map fi — > (rj(fx), h(fJ,), Z2 (/•*)> hit*)) ^ s regular at fx — 0. 

Then, by the introduction of a complex variable, the above system reduced to 
the family of parameter-dependent continuations of the center manifold, is 
orbitally topologically equivalent to 

w' — (77 + iuio)w + tvj\w\ 2 + uw\w\ A + aw\w\ 6 + Uw\w\ 8 

where r\, r, v and a are unfolding parameters. 

Remark 1.4 The expressions for the Lyapunov coefficients in this article 
will be of interest in the study of the Hopf bifurcations of codimensions 1, 
2, 3 and 4 in the Watt governor system with a spring JEj, pursuing previous 
bifurcation analysis carried out in J?J/. 
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